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1. Introduction 

As it was shown in a number of papers (see [1-5] and references herein), the Poincare gauge 
theory of gravity (PGTG) offers opportunities to solve some principal problems of general 
relativity theory (GR) and modern cosmology. It is because the gravitational interaction in 
the framework of PGTG can have the repulsive character under certain physical conditions 
in usual gravitating systems with positive values of energy density and pressure. As a result 
the PGTG allows to build totally regular Big Bang scenario with an accelerating stage of 
cosmological expansion at present epoch. Such scenario is based on homogeneous isotropic 
models (HIM) with two torsion functions built in the framework of PGTG in our previous 
paper [4], where as an example a particular regular inflationary cosmological solution was 
obtained by numerical calculations. 
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The present paper is devoted to analysis of regular inflationary cosmological HIM with 
two torsion functions in the PGTG. In Section 2 cosmological equations for such models in 
dimensionless form are given. In Section 3 the most important general properties of regular 
inflationary cosmological solutions are analyzed numerically. In Section 4 we study extremum 
surfaces, on which the Hubble parameter vanishes and initial conditions for investigated 
solutions are given. 

2. Cosmological equations for HIM with two torsion functions in dimensionless 
form 

In the framework of PGTG any HIM is described by three geometric characteristics: the 
scale factor of Robertson- Walker metrics R and two torsion functions Si and S2 as functions 
of time t. By using general form of gravitational Lagrangian including both a scalar 
curvature and various invariants quadratic in gravitational field strengths — curvature and 
torsion tensors, cosmological equations for HIM with two torsion functions were deduced 
in [4]. The energy density p and pressure p of gravitating matter play the role of sources 
of gravitational field in cosmological equations. These equations contain three indefinite 
parameters: parameter a with inverse dimension of energy density determining the scale of 
extremely high energy densities, parameter b with dimension of parameter / — (167TG) -1 (G 
is Newton's gravitational constant, the light velocity c = 1) and dimensionless parameter e. 
As was shown in [4], by certain restrictions on indefinite parameters cosmological equations 
for HIM lead to observable accelerating cosmological expansion at present epoch. According 
to obtained estimations the parameter b has to be close to f and the parameter e satisfies 
the following conditions 0<£<1. 

In order to analyze inflationary cosmological solutions we will consider HIM filled with 
scalar field <f> with a potential V and usual gravitating matter (values of gravitating matter 
are denoted by means of index "m"). In this case we have: 



To investigate inflationary solutions we transform cosmological equations obtained in [4] 
to dimensionless form by introducing dimensionless units for all variables and parameter b 
entering these equations and denoted by means of ~ by the following way: 
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where dimensionless Hubble parameter H is defined by usual way H = Rr 1 ^. Numerical 
analysis in Section 3 will be made by choosing quadratic potential V = ^m 2 <j) 2 , then in 
accordance with (1) the transition to dimensionless units gives: 

1 

2 

Because of the transformation (1) cosmological equations ((22)-(23) in [4]) take the following 
dimensionless form, where the differentiation with respect to dimensionless time t is denoted 
by means of the prime and the sign of ~ is omitted below: 
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(Z = 1 + p - 3p - 2 (b + e) Si = 1 + AV - <p' 2 + p m - 3p m - 2 (b + e) Si). The torsion function 
Si in dimensionless form entering (2)-(3) is 
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and dimensionless torsion function S'2 satisfies the following differential equation of the second 
order: 

e [S' 2 ' + 3HS' 2 + 3H'S 2 - 4 (S[ - 3HS 1 + AS 2 ) S 2 ] 

-2 (p - 3p - 2bSi) S 2 -2(l-b)S 2 = 0. (5) 

By using dimensionless units the equation for scalar field and conservation law for gravitating 
matter have the usual form 

<f>" + 3H<j>' = ~, (6) 

p' m + 3H(p m + Pm )=0. (7) 

By means of (6)-(7) and (4) for SVfunction we can transform cosmological equations (2)-(3) 
and (5) for ^-function to the following form: 
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The cosmological equation (8) leads to the following equation for extremum surface in space 
of variables (0, 0', S 2 , S' 2 , p m ), in points of which H = 0: 
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where variables on extremum surface are denoted by means of index 0. Then from (9) 
we obtain the following expression for derivative of the Hubble parameter H' in points of 
extremum surfaces 
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The transition from compression stage to expansion stage takes place on extremum surface. 
In the case of HIM filled at the beginning of cosmological expansion with scalar field with 
quadratic potential V and ultrarelativistic matter (p m = |p m ) analyzed below the equation 
(11) for extremum surface takes the following form: 
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3. Numerical analysis of regular inflationary cosmological solutions 

We will obtain cosmological solutions by integrating the system of differential equations (9), 
(10), (6), (7) and by choosing initial conditions for independent physical variables given 
on extremum surface in accordance with (13). Because the most important properties of 
cosmological inflationary solutions at the beginning of cosmological expansion are connected 
with the presence of scalar fields, at first we will analyze HIM filled by scalar field without 
other gravitating matter (p m = 0). For simplicity we will consider flat HIM (k = 0). By 
taking into account restrictions on indefinite parameters leading to cosmological acceleration 
at asymptotics [4], we will use b = 1 for numerical calculations. 

Any regular inflationary cosmological solution includes the following stages: the 
compression stage, the transition from compression to expansion, the inflationary and post- 
inflationary stages. Properties of various stages of regular inflationary cosmological solution, 
generally speaking, depend on parameters e and m and initial conditions for independent 
physical variables. Because of relation (13) only three variables from the following four 
quantities (0o, <f>b, S 2 o, S' 20 ) are independent. We will give the initial conditions at the 
moment t = for (0 O , S 20 , S' 20 ), then the initial value of derivative 0' o will be determined 
from (13). At first, we will study the transition stage from compression to expansion. The 
temporal behaviour of the Hubble parameter H(t), scalar field <fi(t), torsion function S 2 (t) 
and its derivative S' 2 (t) at this stage are presented in Figures 1-3. The graphs in Figures 1-3 
are obtained at the following values of parameters e = 10~ 5 , m — 1.1 and initial conditions: 
0o = 10, 5*20 = and S 20 = 10 -4 . The characteristic feature of transition stage is essentially 
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non-linear oscillating behaviour of the Hubble parameter. There is the correlation between 
the frequency of such oscillations and that of the function S' 2 (see Fig. 1 and Fig. 4). 

H, S 2 




t 



Figure 1. The temporal behaviour 

of H(t) (firm line) and S' 2 (t) Figure 2. The temporal behaviour 

(dotted line) at the transition stage of <j)(t) at the transition stage. 
(S'2 is normalized by the factor 0.1). 




of the transition stage. 



The oscillating behaviour of the Hubble parameter at transition stage allows to estimate 
its duration only approximately. As numerical analysis shows, the duration of transition 
stage is smaller in one order or more than duration of inflationary stage. At the end of 
transition stage the amplitude of if-oscillations decreases with time, the value of H becomes 
positive, and at some moment the transition to inflationary stage with slow rolling regime 
of scalar field takes place. The inflationary stage is finished at t cn d, when the scalar field 
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becomes equal to zero and then it becomes to oscillate. Because, as was noted above, the 
duration of inflationary stage is greater in one order or more than duration of transition stage, 
we can consider the value of t en d as estimation of duration of inflationary stage. Numerical 
investigation of dependence of time t e nd on parameters m, e and on initial conditions for 
0o, 5*20, S20 leads us to the following results. The value of t en d depends essentially only on 
parameter m and initial value of (f> , and by given values of (f> and m it depends weakly on 
admissible values of e and initial values of S 2 q, S' 20 (Figures 5-9). 





Figure 7. The dependence of t en d 
on e (e = 10 £l , (f>o — 10, m = 1.0, 

s 20 = 10- 5 , S' 20 = 0). 



Figure 8. The dependence of tend 
on S'20 (e = 10~ 5 , m = 1.0, 0o = 
25, S' 20 = 10- 4 ). 



Similar to GR, the value of duration of inflationary stage depends on initial value of 
0o in linear way (Fig. 5). From numerical analysis of the value of t en d in dependence on 
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Figure 9. The dependence of i e nd 
on S' 20 (e = 1CT 5 , O = 10, m = 

1.0, s 20 = io- 5 ). 



parameters and initial conditions follows: 

(i) similar to GR, there is a lower limit for initial value of |0o| to have sufficient number of 
e- folds during the inflationary stage (|0 O | ~ 10); 

(ii) for given value of e admissible values for |0 O | have upper limit; 

(iii) the parameter e has upper limit depending on values m, O and S' 20 (in Fig. 7: e < 10~ 3 ); 

(iv) there are upper limits for admissible values of S'20 an d S' 20 at given values of m, e and 



Note that the presence of relativistic matter besides of scalar field leads only to 
quantitative corrections and does not change general conclusions given above. In particular, 
the presence of ultrarelativistic matter at transition stage changes the amplitude of 
oscillations S^-function. The influence of ultrarelativistic matter at inflationary and 
post inflationary stages is negligibly small because the energy density of ultrarelativistic 
matter rapidly decreases during inflationary stage. 

The behaviour of H and during the inflationary stage are similar to that of GR. 
(Figures 10). There are small differences for H at the start and the end of inflationary stage. 
The graphs in Fig. 10 and also in Fig. 11 for ^-function are obtained for m = 1.1, e = 10~ 5 , 
4>q = 10, S 2 o = 0, S' 20 = 10~ 4 . During the inflationary stage the following relations are 
satisfied with a rather high accuracy 



The behaviour of scalar field 0, torsion function S2 and Hubble parameter H during 
the postinflationary stage are presented in Figures 12-14. The graphs in Figures 12, 13 



0o (Fig. 8-9). 
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Figure 10. The behaviour of 
H (firm line) and <j> (dashed line) 
during the inflationary stage. 



Figure 11. The behaviour of 52 
during the inflationary stage. 



and 14a are obtained at the following values of parameters e = 10~ 5 , m = 1.1 and initial 
conditions: 0o — 10, S20 = and S' 20 = 10~ 4 . Similar to GR during the postinflationary 
stage the scalar field <fi oscillates with decreasing amplitude. The frequency of 0-oscillations 
increases by increasing of the parameter m. The oscillations of the Hubble parameter have 
the character of beats (Fig. 14). Two subsequent pulsation of beats are divided by domain 
of oscillations with small amplitude. Similar to inflationary models without torsion function 
S 2 investigated in [3], for sufficiently large values of the parameter m the oscillating Hubble 
parameter changes its sign (in Figures 14a and 14b the value of m is equal to 1.1 and 0.4 
respectively). In this case the cosmological model vibrates during some time interval after 
inflation. The amplitude of //-beats increases by increasing of parameter m. By decreasing 
of e the frequency of oscillation inside one pulsation increases although the frequency of 
beats oneself practically does not change. There is the correlation between 0-oscillations, 
//-beats and the temporal behaviour of ^-function during the postinflationary stage. When 
instantaneous value of is small in absolute value, then ^-function oscillates near zero and 
H decreases being between two subsequent pulsations. When absolute value of is near 
to its maximum, the //-parameter oscillates inside of pulsation and the absolute value of 
^-function increases. 

As result of this Section, we can conclude that physically interesting solutions exist for 
some restrictions on parameters m, e and for sufficiently large domain of initial conditions. 

4. Extremum surface for initial conditions 

It is interesting to investigate more particularly the structure of extremum surface defined 
by (13), on which initial conditions for regular inflationary cosmological solutions are given. 
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Figure 12. The behaviour of < 
during the postinflationary stage. 



Figure 13. The behaviour of 52 
during the postinflationary stage. 
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Figure 14. The behaviour of H during the postinflationary stage. 



Supposing below p = and k = 0, we will have extremum surface in 4-dimensional space of 
variables (<f>o, 0' o , S20, S' 2Q ). By taking into account results concerning the duration of time 
t en d obtained in previous Section, we will consider 3-dimensional subspace P 3 of variables 
(£20, S' 2Q , 0q), then the relation (13) determines 2-dimensional surface P 2 in P 3 depending 
parametrically on variable O an d also on parameters m and e. The 2-dimensional surface 
P 2 has sufficiently complicated structure and includes a closed cover P 20 and also some 
complicated surfaces surrounding the closed cover P 20 . A P ar t of surface P 2 is presented in 
Figure 15. 

As numerical analysis shows, we obtain regular inflationary cosmological solutions by 
choosing initial conditions on closed cover P 20 . The study of geometrical properties of surface 
P 2 allows to understand some important properties of regular inflationary cosmological 
solutions at transition stage from compression to expansion. With this purpose we will 
consider cross-sections of the surface P 2 with planes orthogonal to coordinate axes S20, S' 20 
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and 0' o . As example, in Fig. 16 corresponding cross-section in the plane S' 20 = 10~ 4 is 
presented (m = 1.1, e — 10~ 5 , O = 10). 




Figure 16. Cross-section of 
Figure 15. Part or extrcmum . „ . . 

extremum surface F 2 q °y tne plane 
surface P 2a . a > _m-4 

^20 — iu • 



By increasing of m and the value of variable 0o linear sizes of the closed cover P 2 o 
increase. By decreasing of e linear sizes of P20 increase only in directions of axes S'20 an d S' 20 
and do not change practically in direction of <f)' Q . 

As it was discussed in previous Section, the Hubble parameter oscillates at transition 
stage from compression to expansion. This means that during this stage the quantities H, 
S2, S' 2 , 0, 4>' are changing by intersecting the extremum surface H = and by changing the 
sign of derivative H' . In order to find domains on closed cover P 20 , where the sign of H' is 
positive and negative, we have to consider the intersection of P 2 q with surface H' = defined 
by (12). The sign of H' on P20 is different from different sides of intersection line. The place 
of intersection of P20 with surface H' = depends on the value of 0o- So, in the case 0o — 10 
the domain of intersection corresponds to small values of S' 20 and to values of S'20 near to its 
maximum. By increasing of 0o the domain of intersection moves in direction of decreasing 
values of S 2 o and increasing values of S' 20 . The presence of domains with positive and negative 
values of derivative H' on extremum surface P20 allows to explain the oscillating character 
of the Hubble parameter at transition stage, that is connected essentially with behaviour of 
S 2 -function (see Fig. 3). After the bounce at t — the S 2 -function quickly increases and 
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reaches its maximum. As a result, oscillating function H(t) reaches now extremum surface 
at value of S20, where the derivative H' Q is negative. Together with changing of ^-function 
the process of oscillations of the Hubble parameter is repeated at decreasing values of 0o (see 
Fig. 2) and as result at decreasing values of S^omax- Oscillations of if-function will continue 
even after this function does not reach extremum surface. The amplitude of iJ-oscillations 
will continue to decrease until the inflationary stage will come. 



5. Conclusion 



As follows from our analysis, regular inflationary cosmological solutions for HIM with 
two torsion functions built in the framework of PGTG are realized at certain restrictions 
on parameters of HIM and at large domain of initial conditions on extremum surface 
H = 0. The presence of pseudoscalar torsion function S2 (t) leads to essential changes of 
considered solutions in comparison with inflationary cosmological solutions for HIM without 
^-function analyzed in [3]. At first of all these changes relate to properties of transition 
stage from compression to expansion and post-inflationary stage and are connected with 
oscillating behaviour of the Hubble parameter. Differences of post-inflationary stages for 
cosmological solutions in considered theory in comparison with general relativity theory can 
lead to quantitative differences by transition to radiation-dominated stage, in particular, to 
differences in anisotropy of relic radiation. This means that the building of perturbation 
theory for scalar fields in considered inflationary cosmological HIM is of direct physical 
interest. 
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